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SUMMARY 


This  report  describes  the  development  of  the  equations  of  motion  of  an 
airborne  towed  vehicle  and  the  evaluation  of  its  lateral-directional  dynamic 
stability.  In  the  study  it  Is  assumed  the  towing  cable  is  dynamically 
stable  and  cable  scopes  of  150  to  400  feet  are  deployed. 

Determination  of  che  character  of  the  lateral-directional  motion  of  a 
HTM-11  and  a  modified  HTM-11  towed  vehicles  employing  the  developed  equations 
of  motion  are  discussed. 
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INTRODUCTION 


This  report  discusses  au  evaluation  of  the  lateral-directional  dynamic 
stability  of  an  airborne  towed  vehicle.  In  so  doing,  the  vehicle's  equations 
of  motion  are  developed.  A  simultaneous  solution  of  the  equations  results 
in  a  sixth  order  polynomial  characteristic  equation.  Evaluation  and  inspec¬ 
tion  of  the  roots  of  the  polynomial  determines  the  dynamic  stability  state 
of  the  vehicle. 

The  equations  of  motion  of  the  airborne  towed  vehicle  are  similar  in 
nature  to  those  equations  derived  for  a  fixed  wing  aircraft  in  flight.  The 
exception  in  the  case  of  the  towed  vehicle  is  that  the  equations  herein 
reflect  the  cable  tension  on  the  body  due  to  deployment  of  various  lengths  of 
cable. 

The  derivation  of  the  equations  of  motion  of  a  fixed  wing  aircraft  in 
flight  has  been  published  in  references  (a),  (b)  and  (c).  Therefore,  a  sunmary 
of  the  derivation  of  those  equations  and  a  development  of  the  equations  pecu¬ 
liar  to  the  towed  vehicle  are  discussed  herein. 

The  vehicle  analyzed  is  assumed  to  be  a  lifting  body,  which  has  its  tow 
cable  attachment  point  at  a  location  fox  rd  of  its  center  of  gravity  similar 
to  the  sketch  contained  in  Figure  1.  It  is  also  assumed  that  the  tow  cable 
is  dynamically  stable.  For  a  typical  mission  it  is  assumed  that  the  vehicle 
is  deployed  to  cable  scopes  of  150  to  400  feet. 

Two  applications  of  the  derived  equations  of  motion  to  real  life,  towed 
body  incidents  are  discussed. 

The  effort  described,  herein,  was  achieved  under  the  authority  of  ref¬ 
erence  (d) . 


EQUATIONS  OF  MOTION  WITH  RESPECT 

TO  AXES  FIXED  IN  SPACE 


It  is  assumed  in  the  analysis  that:  (1)  the  vehicle's  airframe  is  a 
rigid  body  (2)  the  earth  is  fixed  in  space  and  (3)  the  mass  of  the  vehicle 
is  constant. 

Mainly  the  equations  of  motion  of  the  vehicle  are  based  on  Newton's  second 
law  of  motion.  This  law  states  that  the  rate  of  change  of  momentum  of  a  body 
is  proportional  to  the  force  applied  to  chat  body.  Also,  the  rate  of  change 
of  the  moment  of  momentum  of  a  body  is  proportional  to  the  torque  applied. 

The  summation  of  the  forces  and  torques  applied  to  the  body  are  referred  to 
a  right  hand  system  of  Cartesian  axes  fixed  in  space  as  ehown  in  Figure  1. 

Mathematically  stated  these  equations  of  motion  are: 
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efYi 

EL 


dC"wV) 

dt 

dt 


(1) 


where  IFy^.  and  V  are,  respectively,  a  summation  of  the  external  forces 

applied  parallel  to  the  Y  axis  and  a  component  of  the  total  velocity  along  the 
same  axis  and  ny  is  the  vehicle's  mass;  where  L  and  N  are  sunmations  of  com¬ 
ponents  of  moments  about  the  Xi  and  Zi  axes,  respectively;  where  hXj  and  hzj 
are  components  of  the  moment  of  momentum  about  the  Xx  and  Zx  axes,  respectively. 

The  components  of  the  moment  of  momentum  are  derived  by  summing  the  moments 
of  the  linear  velocity  vectors  about  each  axis,  multiplying  by  the  differential 
mass  and  integrating  over  the  entire  mass.  The  resulting  components  are: 


*  PIXX  "  ^XY  ' 

RIXZ 

"  RI2Z  ’  PIXZ  " 

qiyz 

(2) 

where  P,  Q  and  R  are  angular  velocities  along  the  Xx»  Yx>  Zj  axes,  respectively, 
and  I  is  the  moment  and  product  of  inertia  about  the  respective  subscripted 
axes. 

EQUATIONS  OF  MOTION  WITH  RESPECT 
TO  EULERIAN  AXES 


If  the  reference  axes  remain  fixed  in  space  and  the  vehicle  rotates,  the 
moments  and  products  of  inertia  will  vary  with  time  and  functions  of  P,  Q  and 
R  will  appear  in  the  equation  (1)  as  derivatives.  In  order  to  circumvent  the 
use  of  these  terms  an  Eulerian  axis  system  is  adopted.  Therefore,  it  was 
assumed  that  the  X,  Y,  Z  axes,  as  seen  in  Figure  1  with  their  origin  at  the 
vehicle's  center  of  gravity,  are  fixed  to  the  vehicle  and  move  with  it. 

A  gennral  expression  of  equations  (1)  referred  to  the  Eulerian  axes 
become : 
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where  uj  x  VT  it  the  centripetal  acceleration  resulting  from  the  vector  cross- 
product  of  angular_valocity,  u>,  and  the  linear  velocity,  Vr;  and  the  vector 
cross-product  w  x  h  is  the  absolute  rate  of  change  of  moment  of  momentum. 

The  expanded  form  of  equations  (3)  using  the  scalar  covaponents  and 
assuming  the  XZ  plane  of  the  vehicle  is  a  plane  of  symmetry  are: 

SFy  -  (V  +  RU  -  PW) 

EL  -  PIxX  -  RIj£2  +  QR(IZZ  -  lyy)  *  pQ*xz 

EN  -  RIZZ  -  PIXZ  +  PQdyy  -  Ixx)  +  QR1XZ  ,  (4) 

LINEARIZING  THE  EQUATIONS  OF  MOTION 

The  complexity  of  the  analysis  can  be  simplified  by  linearizing  the  equa¬ 
tions  of  motion.  The  equations  are  linearized  by  assuming  that  the  motion  of 
the  vehicle  is  the  result  of  disturbing  it  from  an  equilibrium,  steady  flight 
condition  and  further  assuming  the  perturbations  from  steady  flight  are  small. 
The  velocities  can  then  be  written  as: 

U  -  U0  +  u 

V  «  V0  +  v 

w  ■  w0  +  w 

P  -  p0  +  p 

Q  ■  Q0  +  q 

R  -  Rp  +  r  (5) 

Zero  subscripts  indicate  steady  state  velocities;  lower  case  letters  indicate 
perturbed  quantities.  It  is  also  assumed  the  vehicle  is  originally  in  steady 
flight  with  the  wings  level  and  all  components  of  velocity  are  zero  except 
U0  and  WQ.  Then 


Po  “  Qo  *  Ro  “  Vo  “  0  (6) 

Expanding  equations  (4)  by  incorporating  equations  (5)  and  (6)  and  assuming 
that  the  perturbations  are  small  enough  that  the  products  and  squares  of  the 
velocity  changes  are  negligible,  the  equations  of  motion  reduce  to: 


EFy  “ 

m(v  + 

V  -  woP) 

EL  *■ 

pixx  - 

rlXZ 

EN  - 

rTZZ  ' 

Prxz 

(7) 
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DEVELOPMENT  OF  GRAVITY 
FORCE  COMPONENT 


The  gravity  force  contribute*  components  to  i  sunmation  of  external 
forces  applied  to  the  vehicle.  In  this  report  the  external  forces  are 
limited  to  contributions  along  the  vehicle's  Y  axis.  The  gravity  components 
act  along  steady  flight  axes  which  it  is  assumed  has  an  initial  pitch  angle 
displacement,  0O,  with  respect  to  the  gravity  vector.  Resolving  the  gravity 
force  along  the  steady  flight  axes  produces: 

Wox  "  “  ®W8  9in  0o 


W0y  *  ° 

woz  "  “W*  cos  90 


(8) 


where  g  is  the  gravity  acceleration  constant.  It  is  further  assumed  that  the 
axes  are  disturbed  and  displaced  from  the  steady  axes  by  Eulerian  angles 
ip,  0  <p ,  in  that  order  of  rotation. 


Assuming  that  the  disturbed  Eulerian  angles  and  0O  are  small,  and  the  cosine 
of  the  angles  equal  one  and  the  sine  of  the  angles  equal  the  angle,  the  pro¬ 
ducts  of  these  angles  are  approximately  zero  and  can  be  neglected.  There¬ 
fore,  the  gravity  force  component  acting  along  the  disturbed  Eulerian  Y  axis 
is: 


‘Pmyg 


(9) 


DEVELOPMENT  OF  THE  AERODYNAMIC  FORCE 

AND  MOMENTS 


The  aerodynamic  forces  and  tnments  applied  to  the  airborne  towed  vehicle 
are  due  to  the  resistance  f  the  surrounding  atmosphere  to  the  vehicle's 
motion. 

The  expression  of  the  aerodynamic  force  along  the  vehicle's  Y  axis  can 
be  written  as: 

YF  -  C*  1/2  P  VT%  (1( 

where  CY  is  a  dimensionless  coefficient,  p  is  the  atmospheric  density,  VY  is 
the  velocity  of  the  vehicle,  and  Sw  is  the  vehicle's  wing  area. 
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Tha  rolling  and  yawing  moment  about  tha  vahicla'a  X  and  X  axis  can  ba 
expressed  raapactivaly  aa: 

2 

L  -  Ct  1/2  p  VT  Syb 

and 


N 


2 

1/2  p  VT  Syb 


(ID 


whara  and  Cni  ara  dinanalonlasa  coefficients  and  b  la  Cha  vahicla'a  wing 
span. 

Tha  fore aa  and  moments  acting  on  tha  dlaturbad  vahlcla  can  ba  wrlttan 
as  a  function  of  the  paramatara  with  which  aach  la  known  to  vary  by  expanding 
tha  tarma  in  a  Taylor  aarlaa. 


Therefore,  the  force  Yp  and  momenta  L  and  N  can  be  written  aa: 


Yr 

-  Y„  + 
Fo 

3Y-  r 

+ 

3Y_  i 
3r 

+ 

ilv  v 

3vF 

♦fi** 

+ 

E.  p 
3pF 

■f 

If*1 

L 

"  Lo  + 

—  r 

3r 

+ 

3L  . 

?7r 

+ 

2LV 

3  v 

3L  • 

+  ”T*  V 

3v 

+ 

Skp 

3P 

+ 

N 

"  N0  + 

3N 

r 

3r 

+ 

3N  , 

3r 

+ 

3v 

+  ^  V 

3v 

+ 

3P  V 

+ 

If  p 

(12) 


The  zero  subscript  indicates  the  portions  of  the  force  and  moments  which  ara 
acting  during  steady  flight.  The  partial  derivative  terms  express  tha  forces 
and  moments  created  by  tha  disturbed  parameters.  In  this  case  due  to  steady 
flight  assumptions  YpQ  -  Lo  ■  N0  -  0.  It  is  assumed  that  tha  air  flow 
about  the  vehicle  is  quasi-steady,  therefore  all  derivatives  with  respect  to 
the  rate  of  change  of  velocities  are  omitted. 


DEVELOPMENT  OF  THE  CABLE  TENSION 
FORCE  COMPONENT  AND  MOMENTS 


It  is  assumed  in  this  analysis  that  the  tow  cable  discussed  herein  is 
dynamically  stable. 

The  orientation  of  the  cable  tension  vector  relative  to  the  towed  vahicla'a 
axis  system  is  presented  in  Figure  2. 

Cable  tension  contributes  components  to  the  summation  of  applied  external 
forces  and  moments  and  is  assumed  to  act  along  the  longitudinal  axis  of  the 
cable. 
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The  steady  stete  components  of  the  tension  which  have  an  Initial  pitch 
angle  displacement  0O  and  a  lateral  angle  displacement  of  a  ere  determined 
by  resolution  of  the  tension  vector  along  the  steady  state  axes  Xq,  Yd  and 
Zo. 


The  steady  state  components  of  the  cable  tension  are: 

T^  -  T  coe  9C  -  0O  T  sin  6C 

Trt  ■  -  0  T  cos  0 

Uy  J 

T0_  ■  0  T  cos  0  +  T  sin  0  / 

Z  o  c  c  V 

where  angles  0O  and  cr  are  assumed  to  be  small  and  T  is  the  tension  of  the 
cable  exerted  on  the  drogue  during  trimmed  steady  flight.  The  component  of 
cable  tension  acting  along  the  disturbed  Y  Eulerian  axis  is  determined  by 
performing  an  Euler  an^le  transformation  and  employing  the  small  angle 
assumptions .  Therefore 


(14) 


Expanding  FyT  by  substituting  equations  (13)  and  simplifying  gives: 


Fyr  -  ’  4* T  cos  0 c  "  <JT  C03  9C  +  <j>T  sin  Pc 


(15) 


The  cable  attachment  point  (tow  point)  on  the  towed  vehicle  is  remote 
from  the  vehicle's  center  of  gravity  as  shown  in  Figure  1,  therefore  moments 
about  the  X  and  Z  disturbed  Eulerian  axes  are  generated  by  the  cable  tension 
component,  Fy*.,  in  equation  ^15).  The  moments  are: 


Ly 

*  (  -  <I>T 

C08 

9c  - 

OT  cos  9  c 

+ 

4>T 

sin 

6c> 

ZTP 

Nij 

-  (  -i|)T 

COS 

9c  - 

aT  cos  0  c 

+ 

$  T 

sin 

9c> 

XTP 

(16) 

S  U 

MMAIIO 

N 

0  F 

E  X  T  E  R 

N  A 

L 

F  0 

R  C 

E  S 

A 

N  D 

MOMEN 

T 

The  external  forces  applied  to  the  towed  vehicle  are:  gravity  force, 
aerodynamic  forces  and  cable  force. 

Substituting  the  forces  acting  along  the  disturbed  Eulerian  Y  axis  from 
equations  (9),  (12)  and  (14)  into  the  force  equation  of  equations  (7)  and 
rearranging  results  in  the  following: 
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V  + 

uor 

-  woP  - 

<t>8  - 

Yrr  +  Yvv 

+  YpP  - 

ij)T  coe9c 

CTT  CO80c  + 

$T  sin9c 

“w 

“W 

where 

_1 

3Yf  ,  _i 

3r  ^ 

3YF 

3v 

and  -1  £* 
%  3P 

are  replaced 

by  Yr, 

Yv 

and  Yp,  These  terms  are  further  defined 

in  Table  I. 

(17 


Substituting  the  moments  about  the  disturbed  Eulerian  X  and  Z  axis  from 
equations  (12)  and  (16)  into  the  moment  equations  of  equations  (7)  and 
arranging  the  terms,  results  in  the  following  equations: 


re- 


fxz  " 

rxx 

L  v 

V 

+  L  p  +  L  r 
P  r 

Ixz 

Izz 

N  v 

V 

+  N  p  +  N  r 
P  r 

T  cose  2  0  T  cose 


c  "TP"  A  T  3in6( 


XX 


XX 


XX 


zz 


zz 


zz 
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where  •p-  and  ~  are  replaced  by  L  and  N 

ZS  v  x 


XX 


Similar  substitutions  are  made  for  functions  of  p  and  r.  l  ,  l  ,  L  N  N 
and  are  defined  in  Table  I.  v  r'  v'  p 


RELATIONSHIP  BETWEEN  EULER 
ANGLES  AND  ANGULAR  VELOCITIES 


As  mentioned  previously,  the  Eulerian  axes  are  displaced  from  the  steady 
flight  axes  by  the  Eulerian  angles  <j/,  0  and  tf>.  The  rates  of  changes  of  the 
Eulerian  angles  can  be  represented  as  vectors  along  the  axes  about  which  the 
rotations  take  place. 


The  following  lateral-directional  equations  relating  these  terms  are: 


P  +  Q  tan  9  sin  <f>  +  R  tan  9  cos  <J> 


<|>  -  R 


fsini] 

\cosdJ 


(19) 
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Restricting  the  analysis  to  the  lateral -directional  mode  and  assuming 
that  the  perturbations  evolve  from  a  straight  and  level,  steady  flight  con¬ 
dition.  P  *0  ■  R  *  i b  m  <b  *q*0  and  0  *  a  and  are  of  small  magnitude, 

ooooo  oo 

Equations  (19)  reduce  to 

|  -  p  +  r  tan  6^  (20) 


DEVELOPMENT  OF  TOW  CABLE 
ANGULAR  ROTATION  RATE 


The  absolute  linear  velocity  of  the  towpoint  with  respect  to  the  towed 
vehicle's  center  of  gravity  is  expressed  as: 

TP  eg  TP/  (21) 

eg 

where  V  is  the  linear  velocity  vector  of  the  center  of  gravity  and  Vxp/C.g  i3 
eg 

the  linear  velocity  vector  of  the  tow  point  with  respect  to  the  center  of 
gravity.  Equation  (21)  can  also  be  expressed  as: 

V  =  V  +  a)  x  r 

TP  eg  c  (22) 

Recaliing_that  the  components  of  the  vehicle's  angular  velocity  vector  u)  are^ 

P_,  Q  andJR  and  the  components  of  rc  (distance  between  eg  and  towpoint)  are 
yTP  ZTPf  the  crcss  Product  of  equation  (22)  can  be  expanded  by  the  followxAg 

determinant : 


w  x  r 


X _  y 

TP  TP 


Expanding  the  determinant 

u  x  r  =  i  (-QZjp  -  RyTp)  +  j  (PzTp  +  Rx^)  +  k  (Py^  -  Qx^) 

The  component  of  towpoint  velocity  along  the  vehicle's  Y  axis  is: 

L  =  Y  +  pz_  +  rx _ 

TP  eg  r  TP  TP 


where  as  before  Pq  =  R  =  o.  The  inertial  velocity  component  Y  can  be 
expressed  in  term!  of  disturbed  axis  velocity  components  by  using  the  following 
Euler  angle  trams formation: 
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Y 

eg 


U  cos  0  sin  1J1  +  V  (sin  sf  sin  9  sin  $  +  cos  $  cos  <Ji) 

+  W  (cos  $  sin  0  sin  -  sin  ^  cos  ifO 


(26) 


Recalling  that  4,  rp,  0  are  small,  expanding  the  velocity  components,  as  during 

the  linearization,  assinning  V  is  zero  and  neglecting  resulting  small  terms, 

Y  can  be  written  ass 
eg 


Y  *0i)i  +  v-  W4 

eg  o  o 


(27) 


Substituting  Y 


eg 


into  equation 


v25) 


U  i|>  +  v-W<J>  +  pz^  +  rx_ 

o  o  TP  TP 


(28) 


The  vehicle  is  attached  to  the  towing  aircraft  as  shown  in  Figure  2.  The 
tow  cable  of  length,  Z,  generates  an  angle,  a,  in  the  XY  inertial  plane  when 
the  vehicle  rotates  about  ^he  cable  attachment  point  on  the  aircraft.  The 
inertial  linear  velocity,  Y^,  therefore,  also  can  be  expressed  as: 


Y^  =■  a  i  cos  ec 


(29) 


Equating  the  equations  (28)  and  (29)  and  rearranging,  the  expression  for  a  is: 


a 


U_\j>  +  v  -  W  <b  +  pz_  +  rx__ 

oT  o  •  TP  TP 

*COS0 


(30) 


AXIS  SYSTEM 


A  body  axis  system  is  chosen  where  the  X  axis  is  oriented  along  the  lon¬ 
gitudinal  principal  axis  of  the  vehicle.  The  product  of  inertia,  I  ,  used  in 
equations  (31)  therefore  is  equal  to  zero.  XZ 


LATERAL-DIRECTIONAL  EQUATIONS 

OF  MOTION 


A  recapitulation  of  the  rearranged  equations  of  motion  of  the  disturbed 
towed  vehicle  taken  from  equations  (17),  (18),  (20)  and  (30)  is: 

v  -  Yvv  -  (Y  +  Wo>  p  +  <UQ  -  Yc)  r  -  (g  ^  ^.3 a)  * 

row  > 

T  c°s  8c  +  T  cos  0c  a  -  0 

+  ™W~ 


1 

i 

i 
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-Lv  +  p-Lp  L  r  - 
v  *  —  r 


z_T  sin0 
TP  ■  c 

Ixx 


♦  + 


Z  T  0030 
TP  c 

xxx 


*  + 


Z__T  CO80 
TP  c 

xxx 


x  T  sin0  x_T  cos  9 

V  •  V  +  r  "  Nrr  "  - ■  *  +  - - C  *  + 


zz 


zz 


x__T  cos0 
TP  c 


ZZ 


(31) 


p  +  r  tan,  <J>  *  0 

r  ~  <p  •  0 


-  v  - 


w  -  w  +  v  -  v 


£cos8 


+  0  =  0 


CHARACTER  OF  MOTION  OF 
DISTURBED  VEHICLE 


A  simultaneous  solution  of  the  equations  of  motion  results  in  a  sixth 
order  polynomial  characteristic  equation.  The  character  of  the  lateral- 
directional  motion  of  the  towed  vehicle  deployed  at  various  cable  lengths  is 
revealed  by  the  roots  of  the  polynomial.  If  the  roots  are  real  numbers,  the 
motion  is  aperiodic.  This  motion  is  convergent  if  the  roots  are  negative, 
divergent  if  positive.  When  the  roots  form  a  complex  pair,  the  motion  of  the 
vehicle  is  oscillatory.  Damped  oscillatory  motion  prevails  if  the  real  part 
of  the  root  is  negative,  undamped  motion  if  the  real  part  is  positive. 


VALIDITY  OF  METHOD 


In  order  to  establish  credence  in  the  ability  of  the  equations  of  motion 
and  resulting  characteristic  equation  to  determine  the  lateral-directional 
dynamic  stability  of  a  given  vehicle,  two  applications  of  the  method  are 
discussed. 

During  a  typical  deployment  of  an  HTM-11  airborne  towed  vehicle,  a  drogue 
used  to  carry  electronic  equipment,  severe  lateral,  undamped  oscillations 
occurred.  Oscillations  of  the  vehicle  were  severe  enough  to  separate  the 
towing  cable.  The  incident  occurred  at  approximately  210  KIAS  at  an  alti¬ 
tude  of  10,000  feet.  Cable  deployed  at  the  time  of  the  incident  was  approx¬ 
imately  three  to  seven  feet  in  length, 
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The  pertinent  mass  properties  and  physical  characteristics  of  the 
HTM-11  are  noted  in  Table  II.  Static  aerodynamic  coefficients  which  are 
listed  in  Table  III  are  based  on  wind  tunnel  data  determined  at  the  Naval 
Ship  Research  and  Development  Facility.  Dynamic  aerodynamic  coefficients 
were  derived  at  the  Naval  Air  Development  Center.  The  perturbed  motion  in 
the  application  was  assumed  to  evolve  from  the  derived  trim  conditions  of 
Table  IV. 

A  characteristic  equation  was  ultimately  determined  after  applying  the 
values  noted  in  the  tables  to  equations  (31).  Three  sets  of  complex  roots 
were  evaluated  for  cable  lengths  between  1  and  300  feet.  A  root  locus  of 
the  results  is  presented  in  Figure  3.  The  portion  of  the  locus  appearing 
to  the  right  of  the  vertical  imaginary  axis  indicates  that  unstable  vehicle 
flight  existed.  The  locus  also  revealed  that  an  unstable  natural  frequency 
of  7.5  radians  per  second  was  evident  at  a  cable  scr-pe  of  5  feet.  In  com¬ 
parison,  an  unstable  natural  frequency  of  7.8  radians  per  second  at  a  similar 
cable  length  was  measured  from  a  video  tape  recording  of  the  unstable  inci¬ 
dent  of  the  HTM-11. 

Further  investigation  employing  the  above  method  revealed  that  by  ex¬ 
tending  the  vertical  end  plates  of  the  wing  of  the  HTM-11  downward,  the 
lateral  dynamic  stability  of  the  vehicle  would  be  stabilized.  Aerodynamic 
characteristics  of  the  modified  HTM-11,  determined  in  the  wind  tunnel  fit  the 
Naval  Ship  Research  and  Development  Facility  and  derived  at  the  Naval  Air 
Development  Center  are  listed  in  Table  III.  Derived  equilibrium 
trim  conditions  resulting  from  the  modification  are  shown  in  Table  IV.  Using 
the  above  data  in  the  equations  of  motion  and  ultimately  solving  for  the 
characteristic  roots  at  various  cable  lengths,  the  root  locus  plot  presented 
in  Figure  4  was  determined.  Three  sets  of  stable  complex  roots  representing 
the  characteristic  equation  of  the  modified  HTM-11  can  be  seen. 

Subsequent  flight  tests  of  the  modified  HTM-11  supported  the  findings  of 
the  study. 


CONCLUSION 


The  equations  of  motion  discussed  in  this  report  can  be  used  to  determine 
the  lateral-directional  stability  of  a  lifting  vehicle  towed  behind  an  air¬ 
craft.  It  is  assumed  that  the  tow  cable  is  deployed  to  lengths  up  to  400  feet 
and  is  dynamically  stable.  Employing  the  equations,  the  stability  of  the 
towed  vehicle  can  be  evaluated  at  various  cable  lengths. 
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b 


c 


eg 

cm 

d_ 

dt 


LIST  OF  ABBREVIATIONS  AND  SYMBOLS 
wing  span 

rolling  moment  coefficient 
yawing  moment  coefficient 
side  force  coefficient 

rolling  moment  coefficient  change  due  to  variation  in  yawing  velocity 

rolling  moment  coefficient  change  due  to  variation  in  rolling  velocity 

rolling  moment  coefficient  change  due  to  variation  in  sideslip 

yawing  moment  coefficient  change  due  to  variation  in  yawing  velocity 

yawing  moment  coefficient  change  due  to  variation  in  rolling  velocity 

yawing  moment  coefficient  change  due  to  variation  in  sideslip 

side  force  coefficient  change  due  to  variation  in  yawing  velocity 

side  force  coefficient  change  due  to  variation  in  rolling  velocity 

side  force  coefficient  change  due  to  variation  in  sideslip 

mean  aerodynamic  chord 
center  of  gravity 
centimeter 

derivative  with  respect  to  time 


deg  degree 
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F  force 

fc  feet 

2  2 

g  gravitational  constant,  32.2  ft/sec  (980  cm/sec  ) 

h  moment  of  momentum 

I__  mass  moment  of  Inertia  with  respect  to  vehicle's  x  axis 

AA 

I  product  of  mass  moment  of  inertia  with  respect  to  vehicle's  X  and  Y  axes 

AX 

1^2  product  of  mass  moment  of  Inertia  with  respect  to  vehicle's  X  and  Z  axes 

mass  moment  of  inertia  with  respect  to  vehicle's  Z  axis 

1  unit  vector  along  X  axis  of  vehicle 

j  uuit  vector  along  Y  axis  of  vehicle 

jai  imaginary  segment  of  complex  root 

k  init  vector  along  Z  axis  of  vehicle 

kg  kilograms 

l  tow  cable  length 

Ijs  pounds 

mass  of  vehicle 

m  metre 

P  rolling  velocity 

p  perturbed  rolling  velocity 

Q  pitching  velocity 

q  perturbed  pitching,  velocity 

R  yawing  velocity 

rad  radian 

r  perturbed  yawing  velocity 
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wing  area 
seconds 

linear  velocity  component  along  vehicle's  X  axis 

3 

cable  tension 

perturbed  linear  velocity  along  vehicle's  X  axis 
linear  velocity  component  along  vehicle's  Y  axis 
total  linear  velocity 

perturbed  linear  velocity  along  vehicle's  Y  axis 
linear  velocity  component  along  vehicle’s  Z  axis 
perturbed  linear  velocity  along  vehicle's  Z  axis 

distance  between  vehicle’s  tow  point  and  eg  along  vehicle's  X  axis 

distance  between  vehicle's  tow  point  and  eg  along  vehicle's  Y  iixis 

distance  between  vehicle's  tow  point  and  eg  along  vehicle’s  Z  axis 

trim  angle  of  attack 
sides i ip 

E>  lar  angle  displacement  in  pitch 

cable  angle  referenced  to  equilibrium  relative  wind 

equivalent  to  angle  of  attack  during  level  flight 
summation  of  iurms 

cable  rotation  projected  in  X^Y^  plane 
real  segment  of  complex  root 

Euler  angle  displacement  in  roll 
Euler  angle  displacement  in  yaw 
angular  velocity 
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.  derivative  of  variable  with  respect  to  time 

vector  quantity 

£ 

—  partial  derival ive  of  varicble  with  respect  to  time 

Subscripts 
F  for  > 

I  respect  to  inertial  axis 

T  respect  to  trimmed  flight 

Y  respect  to  Y  axis 

0  initial  conditions 

UNITS 

United  States  (U.S.)  customary  units  were  used  in  this  study,  but  in  the 
reporting,  both  metric  and  U.S.  customary  units  appear.  Metric  units  are 
enclosed  by  parSfthesls. 
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TABLE  I 

LATERAL-DIRECTIONAL  STABILITY  DERIVATIVES 


N. 


pV  S  b 
K  T  W 


21 


ZZ 


pvmS„b 
H  T  W 


21 


XX 


pVmS  b 
K  T  W 


41 


ZZ 


c, 

* 


Y 

r 


pvTswb 


4m 


W 


»v/J 


Y 

v 


DV  s 
y  T  W 


2m 


W 


3 


Y 

P 


pVwb 

4mw 
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TABLE  II 

PHYSICAL  AND  MASS  PROPERTIES 

b  4.0  ft  (1.219m)  Sw  5.50  ft2  (0.51»2) 

c  1.375  ft  (0.419m)  Weight  225.  (lbs)  (102. kg) 

Ivv  5.46  slugs-ft2  (7.40  kg-m2)  x  2.916  ft  (0.889m) 

J__  49.26  slugs-ft2  (66.80  kg-m2)  z_  0.812  ft  (0.248m) 

ZZ  Tr 
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TABLE  III 

LATERAL-DIRECTIONAL  AERODYNAMIC  COEFFICIENTS 


HTMrll 

HTM-11 

(MODIFIED) 

C  (pad) 

B 

C  (rad,) 
n8 

C  (rad) 

a 

C  (rad) 

a 

C  ( 
“b 

=r) 

k 

c 

xe 

-0.155 

0.802 

-2.349 

-0.0287 

1.347 

-2.607 

C£  (rac*) 

P 

C£  (rad) 
r 

C  (rad) 
n  '  ' 

r 

C£  (^5d) 
P 

C*  ^ 
r 

rad/ 

c  (£) 

”r 

-0.272 

0.029 

-1.016 

-0.280 

0.070 

-1.399 

TABLE  IV 

TRIM  CONDITIONS 

a  (deg. 5 

0  (deg.) 
c 

T  (lbs.)  a  (deg.) 

9c  (deg.) 

T  ( lbs . ) 

-1.78° 

-41.11 

246.90  -2.20 

-37.80 

250.83 
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Cable  Angle  Sign  Convention: 

Vehicle  below  towing  aircraft  (-) 

Vehicle  to  right  of  towing  aircraft  a  (+) 
(looking  forward) 


FIGURE  2  ORIENTATION  OF  CABLE  TENSION 


FIGURE  3  ROOT  LOCUS  -  HTM-11  CHARACTERISTIC  EQUATION 
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